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ABSTRACT 

We present a comprehensive study of the thermonuclear bursts and millihertz quasi-periodic oscil- 
lations (mHz QPOs) from the neutron star (NS) transient and 11 Hz X-ray pulsar IGR J17480-2446, 
located in the globular cluster Terzan 5. The increase in burst rate that we find during its 2010 
outburst, when persistent luminosity rises from 0.1 to 0.5 times the Eddington limit, is in qualitative 
agreement with thermonuclear burning theory yet opposite to all previous observations of thermonu- 
clear bursts. Thermonuclear bursts gradually evolved into a mHz QPO when the accretion rate 
increased, and vice versa. The mHz QPOs from IGR J17480-2446 resemble those previously observed 
in other accreting NSs, yet they feature lower frequencies (by a factor ~3) and occur when the per- 
sistent luminosity is higher (by a factor 4-25). We find four distinct bursting regimes and a steep 
(close to inverse cubic) decrease of the burst recurrence time with increasing persistent luminosity. 
We compare these findings to nuclear burning models and find evidence for a transition between the 
pure helium and mixed hydrogen/helium ignition regimes when the persistent luminosity was about 
0.3 times the Eddington limit. We also point out important discrepancies between the observed bursts 
and theory, which predicts brighter and less frequent bursts, and suggest that an additional source 
of heat in the NS envelope is required to reconcile the observed and expected burst properties. We 
discuss the impact of NS magnetic field and spin on the expected nuclear burning regimes, in the 
context of this particular pulsar. 

Subject headings: accretion, accretion disks — binaries: close — globular clusters: individual (Terzan 
5) — stars: neutron — X-rays: binaries — X-rays: individual (IGR J17480-2446) 



1. INTRODUCTION 

Matter accreted onto neutron stars (NSs) is piled up 
and compressed, settling towards regions of increasing 
density and temperature. In this process, depending on 
the rate at which accretion proceeds, both stable and 
unstable thermonuclear burning of the accreted H and 
He in to heavier elements are expected (Fuii moto et al.l 
Il981h . The main parameter thought to determine the dif- 
ferent burning regimes is the mass accretion rate on the 
NS per unit s urface area, m fe.g. iFuiimoto et all 119811 : 
iBildsterJ 119981) . When the burning layer becomes ther- 
mally unstable heat cannot be transported as fast as it 
is produced and a thermonuclear runaway occurs, pro- 
ducing a shell flash that releases lO'^^-lO"^^ erg in tens 
of seconds (we do not discuss herein long bursts and su- 
perbursts, which are mo re energetic, 10^°-10"^^ e rg, and 
much less common; e.g., iKeek fc in't Zand! 120081 ). Most 
of the energy is deposited in the outermost layers of 
the NS within a few seconds and radiated away ther- 
mally for tens of seconds while the photosphere cools 
down. Type I X-ray bursts, which feature the spec- 
tral imprint of such photospheric cooling, were discov- 
ered in low-mass X-ray binaries (LMXBs) more than 
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30 v ears ago (iGrindlav et all 119761: iBelian et all 119761: 
Hoffman et al. 1 11978D and promptly i dentified as ther- 
monuclear bursts from accre ting NSs (iWo oslev fc TaamI 
19761: fMaraschi fc Cavalierd 119771 : iLewin fc Joss 1977|). 
Joss fc Lil ()1980f ) pointed out that a strong NS magnetic 
field can act to stabilize nuclear burning in different ways, 
which may explain the fact that no thermonuclear bursts 
have been observed to date from accreting NSs in high- 
mass X-ray binaries (HMXBs). 

Although direct observational evidence of stable ther- 
monuclear burning on accreting NSs has been elu- 
sive, as this is outshined by the much more ef- 
ficient accretion-powered "persistent" emission, theo- 
retical arguments and indirect obs e rvational evidence 
jTaam" '1981': 'Fuiimoto et al.' '1981'; 'van Paradii s et al.l 
Il988: Xewin et al.. .1993. : .Bildsten ,1998.) suggest that 
at very high ifi, close to or above the Eddington 
mass accretion rate (Sees. [5] & S]), thermonuclear burn- 
ing of the accreted H and He proceeds only stably. 
Near the transition between unstable and stable burn- 
ing, an os cillato ry burning regime was predicted by 
Paczvnski (119830. known as marginally stable burn- 
ing. Revnivtsev et al.l ()2001| ) discovered millihertz quasi- 
periodic oscillations (mHz QPOs) in the X-ray flux of 
three atoll sources (the sub-class of low luminosity NS- 
LMXBs; JHasinger fc van der Klii Il989l): 4U 1636-536, 
4U 1608-52 and Aql X-1 (see also iStrohmaver fc Smith] 
1201 It ). They attributed this new phenomenon to 
marginally stable burning on the NS surface. The mHz 
QPO frequency in one of these systems has been found 
to decr ease with time until a bright type I X-ray burst 
occurs (jAltamirano et aLll2008[ ). The persistent luminos- 
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Fig. 1. — Overview of the burst and persistent emission evolution along the outburst of T5X2. Black lines show total 2 s time resolution 
light curves during one RXTE orbit for nine selected dates, as indicated (using PCU2 and the full ~2-60 keV band; gray dashed line 
shows the approximate, not subtracted, background rate). Times have been shifted arbitrarily for display purposes. The increase in burst 
rate and decrease in burst brightness as the persistent flux rises is evident. The step in October 13 was produced by a lunar eclipse 
llStrohmaver et ani2010l) . The source became unobservable for RXTE after November 19 due to Solar constraints. 



ity at which such mHz QPOs are observed has remained 
a puzzle, as it suggests a critical m about an order of 
magnitude low er than the stability boundary predicted 
by theory (e.g.. lHeger et al.ll2007bl ). 

On 2010 October 10, an X-ray transient in the di- 
rection of the globular cluster Terzan 5 was discov- 
ered with t he International Gamma-ray Astrophysic s 
Laboratory ("Bord as et all 120101: iChenevez et alj I2O1O0 . 
During the following week, Rossi X-ray Timing Ex- 
plorer (RXTE) observations revealed 11 Hz pulsa- 
tions (jStrohmaver fc Markwardtl l2010|) and burst os- 
cillations at the same fre quency ijAltamirano et aP 
I2010bt iCavecchi et alJ [201 ID . The Chandra localiza- 
tion (jPoolev et al.ll201C ) confirmed that this was a new 
NS transient, named ICR J17480-2446 (labeled CX2 5 
or CXOGlb J174804.8-244648 by .Heinke et all [2006I ). 
We refer hereinafter to ICR J 17480-2446 as T5X2, as 
this is the second bright X-ray source discovered in 
Terzan 5 (after EXO 1745-248). A 21.3 hr orbital pe- 
riod was measured from the Doppler shift s on the pulsa r 
frequency (Strohmaycr et al. 2010; Papi tto et al.|[20TlT ). 
and the NS magnetic field was estimated to be between 
10^- 10^° G based on the inferred magne tospheric ra- 
dius (jPapitto et al.ll20 lit [Miller et al.ll2011| ). This makes 
T5X2 the type I X-ray burst source (burster) with the 
slowest known NS spin and with the highest estimates 
of the NS magnetic field strength. Near the outburst 
peak T5X2 showed X-ray spectral and variability be- 
havior typical of Z sources (the sub-class of high lumi- 
nosity NS-LMXBs: lHasinger fc van der K"lislll989( ). when 
it was accreting at abo ut half of the Eddington rate 
dAhamirano et al."2010c'). 

Linares et al. (2010a) argued that all the X-ray bursts 
from T5X2 had a thermonuclear origin, based on the 
evolution of the burst rate. Given the lack of spec- 
tral softening along the tail of many of the bursts 



and th eir short recurrence times, I Galloway fc in't Zandl 
(|2010D suggested that some of the T5X2 bursts were 
type II instead of type I (i.e., accretion- instead of 
nuclear-powered). However, the persistent-to-burst en- 
ergy ratio throughout the October-November outburst 
of T5X2 was typical of type I X-ray bursts, i.e., fully con- 
siste nt with the accretion-to-thermonuclear efficiency ra- 
tio (iL inares et al.' 2011"; iChakrabortv fc B hattacharvval 
l2011ir iMotta et al. 2011]). Furthermore, Linares et al] 
()2011|) measured a smooth evolution of the burst lumi- 
nosity and spectral profiles and put forward a scenario to 
explain the lack of cooling in the faintest bursts, conclu- 
sively identifying all X-ray bursts detected from T5X2 as 
thermonuclear. 

We present a thorough analysis of the mHz QPOs from 
T5X2, inclu ding but not li mited to the ones originally 
reported bv iLinares et al.l ()2010al ). We study the mHz 
QPO frequency evolution and energy- dependent ampli- 
tude, as well as all X-ray bursts from T5X2 detected with 
RXTE while the persistent (accretion) luminosity varied 
along the outburst. Unlike previous studies (Motta et al] 
120111 : IChakrabortv fc Bhattacharvval 1201 It) , we analyze 
the complete sample of RXTE bursts and compare their 
properties to theoretical models of thermonuclear burn- 
ing, along the full range in persistent luminosity (~10- 
50% of the Eddington luminosity). Section [2] gives the 
details of the data analysis, and Section [3] presents the 
main observational results: a smooth evolution between 
bursts and mHz QPOs (Figure [T]), the mHz QPO prop- 
erties in detail and four different bursting regimes during 
the October-November 2010 outburst of T5X2. In Sec- 
tion|4]we place the unique mHz QPO and bursting behav- 
ior of T5X2 in the framework of thermonuclear burning 
theory and discuss the possible effects of composition, NS 
spin and magnetic field on the observed bursting prop- 
erties. Section [5] gives our summary and conclusions. 
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2. DATA ANALYSIS 

We analyzed all RXTE observations of T5X2 dur- 
ing its October-December 2010 outburst: a total of 
46 observations taken between 2010 October 13 and 
2010 November 19 (proposal-target number 95437-01). 
The source became Sun constrained after that date, 
and was not detected by the Monitor of All-sky X-ray 
Image (MAXI) on 2010 December 28, indicating that 
the outburst finished between 2010 November 19 and 
2010 December 28. We visually searched for X-ray 
bursts in the full dataset, using 2 s time resolution 2- 
30 keV lightcurves. We performed time-resolved spec- 
troscopy of all bursts using high time resolution data 
(E_125us_64M_0_ls, or GoodXenon when available). We 
extracted dead-time corrected spectra in 2 s time bins, 
using a ~100 s-long pre- or post-burst interval as back- 
ground. We added a 1% systematic error to all channels, 
grouped them to a minimum of 15 counts per channel 
when necessary and fitted the resulting spectra within 
Xspec (v 12.6.0q), using a simple blackbody model with 
the absorbing column density fixed to 1.2x10^^ cm~^ 
(jHeinke et all 120061. We used a distance to T5X2 of 
6.3 kpc, the high est value reported fro m HST photom- 
etry of Terzan 5 (jOrtolani et al.l 120071 ). consistent with 
the distance measurement based on photospheric radius 
expansion bu rsts from another bur ster in the same glob- 
ular cluster (jGallowav et al.|[2008[ ). We note, however, 
that recent estimates of this distance range between 
4.6 kpc and 8.7 kpc (iCohn et al.l 120021 : lOrtolani et all 
120071 : iLanzoni et al.l 120101 ) . and therefore any value of 
the luminosity, energy and mass accretion rate has a sys- 
tematic uncertainty of a factor '^3.6. 

We measured the burst rise time {t-dsc, as defined in 
iGallowav et al.l [20081 i.e., the time to go from 25% to 
90% of the peak count rate), and the total radiated en- 
ergy (i?b) by integrating the bolometric luminosity along 
each burst. We defined the wait time, iwait, as the 
time elapsed between the peak of the previous burst and 
the peak of a given burst, available when no data gaps 
were present before the burst. We measured the daily- 
averaged burst recurrence time, tree, as the total expo- 
sure time during one day (excluding those orbits where 
no bursts are detected) divided by the number of bursts 
detected on that day. Therefore, the instantaneous and 
daily-averaged burst rate, t'burst, are given by ^"2;^ and 
respectively. When only one burst was detected on a 
given day, we considered <ioc an approximate lower limit 
on the recurrence time. We also obtained daily averages 
of -Eb, peak burst luminosity (ipcak), blackbody tem- 
perature (kTpeak) and radius (i^peak )i follow ing the same 
method described in iLinares et al.l (|20111 ): peak burst 
values correspond to a 4 s long interval around the burst 
peak. 

Due to the smooth evolution from a series of bursts into 
a mHz QPO and vice versa (see Fig. [1] and Sec. [3]) the 
distinction between "frequent bursts" and mHz QPOs is, 
to some extent, arbitrary. We searched for mHz QPOs 
all observations taken when the daily-averaged burst re- 
currence time was shorter than 350 s, which corresponds 
to 10 observations between MJDs 55485-55490 (around 
the outburst peak) . Given the typical duration of a con- 
tinuous RXTE observation segment (an "orbit"), this 
threshold ensures that about 10 or more QPO cycles. 



or bursts, are observed without interruption. To study 
the mHz QPOs in those RXTE orbits we used 2-60 keV 
1 s-bin light curves from all active PCUs combined. We 
then calculated a Lomb-Scar gle per i odogram (LSP, over - 
sampled by a factor of 3; ILombl 119761: IScargld [1981 
for each light curve and measured the mHz QPO fre- 
quency, J^QPO) as that frequency with the highest power 
in the periodogram. The corresponding period was then 
used to fold the background-corrected light curve (us- 
ing background rates from pcabackest in 16-s steps) and 
produce a mHz QPO folded profile, from which we mea- 
sured the fractional root-mean-squared (rms) amplitude. 
To investigate the energy dependence of the mHz QPO 
amplitude, we also produced light curves in five differ- 
ent energy bands (in keV: 2-3, 3-5.5, 5.5-9.5, 9.5-21 
and 21-53), folded the light curve at the period found in 
the respective 2-60 keV range dataset and measured the 
"rms spectrum" of the mHz QPOs. We also performed 
2048 s-long fast Fourier transforms (FFTs) using the 2- 
60 keV energy band and the same RXTE orbits, in order 
to constrain the mHz QPO coherence or "quality factor" . 
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Fig. 2. — From top to bottom panels, evolution along the out- 
burst of persistent 2-50 keV luminosity; ii) burst rate as mea- 
sured from (wait (open gray circles), tree (filled black circles) and 
mHz QPO frequency (red filled squares; see Sec. [2] for definitions) 
and in) burst rise time and duration. Gray and black circles show 
individual burst measurements and daily averages, respectively. 
Open triangles show burst rate daily averages based on one sin- 
gle burst, which we consider as upper limits. 

In order to measure the persistent (accretion) luminos- 
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TABLE 1 

Daily- AVER AGED burst properties and persistent luminosity from T5X2. 
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± 


0, 


,13 


4.00 


± 0.28 


55518 


3407 


2 


1703.5 


1, 


,22 


± 


0.03 


3.80 


± 


0.26 


2, 


.31 


± 


0.04 


3, 


,24 


± 


0, 


,11 


4.00 


± 0.30 


55519 


6788 


2 


3394.0 


1, 


,45 


± 


0.04 


3.90 


± 


0.34 


2, 


.19 


± 


0.05 


3, 


,61 


± 


0, 


,16 


3.80 


± 0.27 



''MJD 55482 is 2010 October 13, and MJD 55519 is 2010 November 19. 

''Total daily exposure time in data segments (orbits) where bursts are detected. The total daily exposure time for segments without 
detected bursts is indicated between square brackets (whenever this time is larger than the burst recurrence time). Thus square brackets 
indicate periods of intrinsic burst cessation. 

'^Total number of bursts detected per day. Square brackets indicate periods of intrinsic burst cessation, when no bursts were detected 
despite long enough exposure time. 

•^Daily-averaged burst recurrence time. 

•^Bolometric integrated burst energy (i?b) ^'nd peak burst luminosity (-Lpeak)- Persistent luminosity (L2— 5o) in the 2-50 keV energy band 
(see Sec.[2]for bolometric correction). Peak burst blackbody radius (i?peak) not color- or redshift-corrected. All use a distance of 6.3 kpc. 
^Peak burst blackbody temperature. 



ity, we extracted one dead-time-corrected spectrum per 
observation from Standard 2 data, excluding all bursts 
and subtracting the background spectrum estimated 
with the bright source background model and pcaback- 
est (v. 3.8). We then fitted each persistent spectrum 
with a model consisting of a disk blackbody, a power law 
and a Gaussian line with energy fixed at 6.5 keV, correct- 
ing for absorption as above. We calculated the 2-50 keV 
persistent luminosity (L2-50) from the best fit model. 
Furthermore, we measured the 0.01-2 keV unabsorbed 
flux from a simultaneous fit to the T5X2 spectrum mea- 
sured by Swift-XRT (0.5-10 keV) and RXTE-FCA (2.5- 
25 keV) on MJD 55501 (2010 November 01), extrapolat- 
ing the phabs*simpl(bbody + diskbb) best fit model down 
to 0.01 kcV, and found a bolometric correction factor 



of 1.13. This bolometric correction factor converts 2- 
50 keV into 0.01-50 keV flux, which we take as bolomet- 
ric flux given that the persistent spectrum remains soft 
(photon index 2.4-3.3) throughout the outburst and we 
therefore do not expect sizeable emission above 50 keV 
(with the only exception of the first observation, when 
T5X2 was in the hard state and the photon index was 
'--^1.7). Using the same procedure, we found a bolometric 
correction factor of 1.02 for the 0.5-50 keV lu minosity 
from Cir X-1 reported bv ILinares et al.l |2010b, used in 
Sec. 2]). We applied these bolometric correction factors 
in order to estimate the bolometric persistent luminosity: 
Lpors- We use throughout this work an Eddington lumi- 
nosity of LEdd=2.5 X 10'^® erg s~^ to calculate Eddington- 



normalized L 



pcrs ■ 



To convert from 



pers 



to m, we as- 



5 



sume homogeneous accretion onto the NS and that the 
(general-relativistic) gravitational energy of the in-fahing 
matter is fully converted into radiation at the surface of 
a 1.4 Mq mass, 10 km radius NS. We define the Edding- 
ton mass accretion rate per unit area, mEdd, as the mass 
accretion rate needed to sustain a luminosity equal to 
-^Edd- With this definition mEdd=l-2 x 10^ g cm^^ s^^. 

3. RESULTS 

We found and analyzed a total of 398 X-ray bursts oc- 
curring between 2010 October 13 and 2010 November 19, 
including the faint and frequent bursts near the peak of 
the outburst that form the mHz QPOs (see below). We 
show in Figure [1] an overview of the joint evolution of 
persistent emission and burst properties along the T5X2 
outburst. Bursts become more frequent and fainter when 
the persistent luminosity, Lpors, increases, and they turn 
into brighter and more frequent bursts during the out- 
burst decay. Most interestingly, the X-ray bursts grad- 
ually develop into the observed mHz QPO during the 
outburst rise, and the mHz QPO mutates into a series 
of bursts along the outburst decay (Figure [Ij, an un- 
precedented phenomenon among thermonuclear bursters. 
This also makes T5X2 the most prolific source of ther- 
monuclear bursts known to date (an average burst rate 
over more than a month of 4.9 hr^^; 5.5 hr~^ excluding 
orbits where no bursts were detected), with the short- 
est recurrence times between therm onuclear bursts ob- 
served to date (as short as ~200 s; c.f. lLinares et al.ll2009l : 
IKeek et al.llMol) . 

The smooth metamorphosis between bursts and mHz 
QPOs can be seen qualitatively in Figure [1] and quanti- 
tatively by studying the evolution of several burst prop- 
erties (Figure [2] & Tabled]). Burst rate, rise time and 
duration as well as peak burst lu minosity, Lppak, and to- 
tal radiated energy E^; (see also lLinares et al 11201 It) , all 
evolve gradually while Lpcrs changes by a factor of ~5 
along the outburst. Figure [2] (middle panel) also shows 
that the burst rate equals the mHz QPO frequency dur- 
ing the outburst peak, i'burst='^QPOi as expected given 
that the mHz QPOs are simply formed by a series of 
faint and frequent bursts. This unique behavior makes 
the distinction between bursts and niHz QPOs somewhat 
arbitrary. As explained in Section [2l our practical defini- 
tion of niHz QPO requires burst recurrence times shorter 
than 350 s, so that typically 10 or more QPO cycles are 
observed without interruption. In a few occasions we find 

TABLE 2 

Properties of the mHz QPOs from T5X2. 





OBSID'* 


"QPO 




mBQpo 






/lO* 


(MJD) 




(mHz) 




(%) 


(org/s) 


(g cm 




55485.46 


04-00 [1] 


2.75 ± 0.2 


1 


9 ± 0.1 


9.1 ± 0.4 


4.5 


± 0.2 


66485.63 


04-01 [1,2,3] 


3 ± 0.5 


2 


2 ± 0.1 


9.9 ± 0.4 


4.9 


± 0.2 


55487.43 


06-000 [1,4] 


4.2 ± 0.2 


1 


3 ± 0.1 


11.4 ± 0.4 


5.6 


± 0.2 


55487.62 


06-00 [1,2] 


4.2 ± 0.2 


1 


7 ± 0.2 


12.1 ± 0.5 


6.0 


± 0.2 


66488.26 


07-00 [1] 


4.0 ± 0.5 


1 


4 ± 0.1 


12.0 ± 0.5 


5.9 


± 0.2 


66489.59 


08-00 [1] 


3.75 ± 0.2 


2 


1 ± 0.1 


9.6 ± 0.5 


4.8 


± 0.2 


55490.63 


09-00 [1,2] 


2.9 ± 0.2 


2 


2 ± 0.1 


9.1 ± 0.4 


4.5 


± 0.2 



^For comparison, the reported values for the "low-Lpors mHz 
QPOs" are in the following ranges: !.'qpo=[7-14.3] mHz; 
rmsQpo =[0.7-1.9]% (2-5 keV); Lpers=[0.6-3.5] X lO^'^erg/s 
(Sec. 1221 IRevnivtsev et al.|[200ll: lAltamirano et al.ll20n8r) . 

'QObservation (from proposal-target 95437-01) where the mHz 
QPOs are detected, orbit numbers used indicated between square 
brackets. 



that one burst is missing from the series of regular bursts, 
and the corresponding values of :^burst = ^wait individ- 
ual bursts are a factor ^2 lower than the general trend, 
as can be seen in Figure [2] (middle panel) . Interestingly, 
a similar behavior with sporadic "missing b ursts" is seen 
in the mHz QPO simulations presented in iHeger et al.l 
(|2007bD . 



Oct 18 — 
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Time (s) 

Fig. 3. — Light curves of the mHz QPOs on October 16 (gray) 
and 18 (black), in the 2—60 keV energy band, using 1-s time bins 
and no background subtraction. One RXTE orbit is shown in each 
case, and time is relative to the start of that orbit (2010-10-16 
14:33:27 UTC and 2010-10-18 15:09:55 UTC). The persistent lumi- 
nosity was higher on October 18 (Lpcrs— 0.45-LEdd) than October 
16 (Lpcrs— 0.38LEdd)- To illustrate this, count rate is normalized 
to the number of active PCUs. 




Frequency (mHz) 

Fig. 4. — Lomb-Scargle periodograms of the mHz QPOs on Octo- 
ber 16 (gray) and 18 (black), showing the QPO harmonic structure 
and the change in uqpq. 



3.1. mHz QPOs 

We report the discovery of several instances of mHz 
QPOs during the peak of the T5X2 outburst, on 2010 
Octo ber 16, 18, 19 , 20 and 21 (MJDs 55485-55490; 
Linar es~et al.ll20l"0al for the initial report of mHz QPOs 
on 2010 October 18 and 19). Two examples of mHz QPO 
light curves are shown in Figure |3l each spanning one 
RXTE orbit. Table [2] shows the main QPO properties: 
fractional rms amplitudes between 1.3% and 2.2% (in the 
2-60 keV band) and ^qpo in the range 2.8-4.2 mHz (note 
that the lower end of this frequency range corresponds to 
our mHz QPO definition. Sees. [5] & [3]). Figure U presents 
LSPs of two cases, on October 16 and 18, clearly showing 
the change in t'qpo as well as the harmonic structure (up 
to 4 overtones are visible in the LSPs) . By inspecting the 
power spectra obtained from 2048-s-long FFTs, we find 
that the mHz QPO power is in all cases spread over one 
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- 2600 



2550 
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2450 a 
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Fig. 5. — Folded background-subtracted 2-60 keV light curves of 
the mHz QPOs on October 16 (gray) and 18 (black). Peaks are 
at phase by definition. Two QPO cycles are shown for clarity. 
Two PCUs were active on October 16, yielding a higher collected 
rate. The dashed horizontal line shows the average count rate, 
and the vertical range corresponds to 10% of that value in each 
case, showing that the fractional amplitude was higher on October 
16 (Table [2]|. The net peak burst luminosity was in both cases 



~1. 3x10^7 erg s"! (Table[TJ. 
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Fig. 6. — "RMS spectrum" of the mHz QPOs on October 16 
(gray) and 18 (black). Their fractional rms amplitude increases 
with energy up to ~20 keV. No mHz QPOs are detected above 
that energy, empty triangles show upper limits on their fractional 
rms amplitude. 



or two frequency bins, from which we derive a hmit on 
the fuU-width-at-half-maximum, FWHM<1 mHz. For 
the 2.8-4.2 mHz QPO frequencies this corresponds to a 
lower hmit on the quahty factor (Q=i^qpo/FWHM) of 
Q>3, which reveals a fairly coherent QPO (a fairly con- 
stant twait)- 

We present in Figure [5] light curves folded at the 
mHz QPO period, showing folded burst/QPO profiles 
for the same two dates, October 16 and 18. In both 
cases the burst/QPO profile is peaked and highly non- 
sinusoidal (as seen also in the raw, unfolded light curves; 
Fig. ISj, which explains the harmonic content. On Oc- 
tober 18, when z^qpo was highest (Table [5]), we find a 
nearly symmetric burst/QPO profile. On October 16 
the burst/QPO profile is slightly asymmetric, with the 
rise faster than the decay. We measure a mHz QPO frac- 
tional rms amplitude in the range 1.3-2.2 %, in the total 
(~2-60 keV) PCA band. Furthermore, from our mea- 
surements of the mHz QPO amplitude at different ener- 
gies (Sec. m Fig. ini) we find that its fractional rms am- 
plitude increases between ~2 and ~20 keV, from ~1.4% 
to ~2.8% in the two cases presented in Figure [HI At en- 
ergies higher than 20 keV the mHz QPO is not visible 
by naked eye in the light curves, nor it is conclusively 
detected using LSPs, FFTs or folded light curves. The 



presence of red noise with variable strength gives rise to 
different upper limits on the fractional rms amplitude 
of the mHz QPO above 20 keV (between 1.4% and 3%; 
Fig.[S]). As discussed in Section H?^ such increase in QPO 
amplitude up to 20 keV is in contrast with mHz QPOs 
from other bursting NSs, which showed a fractional am- 
plitude decreasing with increa, sing energy between 2 and 
5 kcV HRevnivtsev et al.ll200lD . 
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Fig. 7. — Burst peak luminosity (top), integrated burst energy 
(middle) and burst recurrence time, tree (bottom) vs. bolometric 
persistent luminosity (Lpcrs). Gray and black circles show indi- 
vidual burst measurements and daily averages, respectively. Red 
filled squares show recurrence times measured from the inverse of 
the mHz QPO frequency. The top axis shows the inferred m as- 
suming homogeneous accretion onto a 1.4 Mq, 10 km radius NS 
(Sec.0. 



TABLE 3 

Broken power law fits to the burst properties vs. 
persistent luminosity relations. 



Relation'' 


-^pcak"-^pcrs 




ircc~-^^pcrs 


[range] 


[0.1-0.3] 


[0.1-0.3] 


[0.17-0.45] 


Ki 


2.2 ± 0.7 


0.9 ± 0.1 


10 ± 7 


ii 


-0.3 ± 0.2 


-0.2 ± 0.1 


-3.2 ± 0.5 


-^'brcak{-f'Edd) 


0.20 ± 0.01 


0.20 ± 0.01 


0.33 ± 0.04 


12 


-2.7 ± 0.1 


-4.4 ± 0.1 


-1.0 ± 0.9 


/ dof 


38.5 / 25 


930 / 25 


6.1 / 30 



''Pairs of variables fitted with a broken power law: Lpcaki 
and iroc 

VS. -L/pcrs 

(see Figures [7] & [Sll . The Lpers range used in 
each fit is quoted between brackets, in units of X/Bdd- ii a-iid 12 
represent, respectively, the power law indices before and after the 
break (at Lbrcak)- 
''Power law normalization before the break, same units as Lpcak 
(lO^'^ erg s-i), Bb (10=*^ erg) and t^oc (s). 



3.2. Burst properties vs. accretion rate 

We present in this Section the relation between burst 
properties and to, as inferred from Lpers, including the 
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Fig. 8. — Same as FigureO showing also the broken power law fits 
to the I/peak I ^bi irec VS. -Lpers relations (Table [Sj. Three regimes 
are apparent and labeled along the bottom axis: slow decrease (at 
ipers/-LEdd$0.2; regime A), fast drop (at 0.2<Lpors/-LEdd$0.3; 
regime B) and saturation (ipers/iEddS^O.S; regime C) of the 
^peaki VS. -Lpers relations. The tree values follow approxi- 
mate tree CX LipcTS 

and tree cx Ljpcrs relations in regimes B and C, 
as indicated (see Figure [TT] Table g] and Sees. I3.2l fcl4t. 
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• At 



the lowest 



persistent 



luminosities, 

0.1<Lpcrs/iEdd<0.2, when Lpers increases both 
Lpcak and i?b (and possibly t^cc) decrease mod- 
erately, from Lpeak— 4.6xl0^''erg/s, i?b— 1-6x10^^ 
erg to ipcak^3.7xl037erg/s, £;b~1.3x lO^^erg. We 
refer to this 0.1-0.2 LEdd regime as slow decrease, 
or regime A (see Fig. [7]). 



• At 



higher 



persistent 



Fig. 9. — Daily-averaged accretion-to-burst energy ratio {a, top 
panel) and burst time scales (rise time and duration, bottom panel) 
vs. persistent luminosity (Lpors) during the T5X2 outburst. Open 
triangles represent lower limits (one single burst detected during 
that day). The corresponding bursting regimes are indicated with 
dotted lines and labels along the bottom axis (Table |4]l. 

faintest and most frequent bursts which form the mHz 
QPOs (see above). Figure [2] shows Lpoak, -Eb and tree 
as a function of the Eddington-normalized Lpors and m 
(Sec. [2]). The overall trend is that of an anti-correlation 
between Lpeak, E^, and tree on the one hand and Lpsrs 
on the other (se e also Fig. 1 in iLinares et al.l 120111 : 
iMotta et al.ll2011[ ). Closer inspection of the burst proper- 
ties over the full Lpers range (0.1-0.5 LEdd; see Figure[7l) 
reveals a more complex and interesting behavior, namely, increases from ^60 to ^120 with increasing Lpers until 



luminosities, 

0.2<Lpers/LEdd<0.3, we find Lpeak and L'b 
to be steeply anticorrelated with Lpers, while tree 
also decreases from 200 s to 400 s with increasing 
Lpers- We refer to this 0.2-0.3 Lsdd regime as fast 
drop, or regime B (see Fig. [T]). 

• At the highest persistent luminosities 
(Lpers>0.3 Lsdd), whcu mHz QPOs are detected, 
burst peak luminosity and energy reach approxi- 
mately constant values {saturation, or regime C): 
Lpeak^l-2xl037erg/s and L;bCi0.3x lO^^erg. The 
burst recurrence time, however, keeps decreasing 
with increasing Lpers from tree— 400 s down to 
-240 s on 2010 October 18 (Fig. (T] Table[2]). 

• Finally, during incursions into the flaring/normal 
branches near the outburst p eak, when T5X2 
showe d Z source behavior (|Altamirano et aH 
l2010cD . neither bursts nor mHz QPOs were de- 
tected. These few episodes without bursting ac- 
tivity, which we refer to as regime D, were as- 
sociated with moderate (~30%) short (hour-long) 
drops in Lpers and occurred between MJDs 55486 
and 55496 (while Lpcrs^O.3 LEdd)- They are noted 
with square brackets in Table [1] 

Interestingly, regimes C and D show similar Lpers (0.3- 
0.5 LEdd), but their variability and spectral properties 
differ (i.e., they constitute diiTerent "accretion states"; 
see Altamirano et al. 2012, in prep.). We summarize the 
main bursting properties of all four regimes in Table |H 

In order to quantify the anticorrelations explained 
above and to characterize the bursting regimes present 
in T5X2, we fit the Lpeak, L^b, tree vs. Lpers rela- 
tions with broken power law functions. The results are 
shown in Figure |S] and Table [31 and confirm that the 
two "breaks" or transitions occur at Lpers— 0.2 LEdd and 
Lpers— 0.3 LEdd- Duc to thc large scatter some of the 
fits are statistically poor, yet they allow us to constrain 
the slope and transition luminosity of the three bursting 
regimes (Fig. El). Remarkably, the ircc-Lpers relation that 
we find in regime B (Lpers~0.2-0.3 Lsdd) is far from lin- 
ear, as would be ex pecte d from a tree cx m^^ relation if 
m (X Lpers (See Sec. l4.1l for further discussion). Instead, 
it is close to tree oc L~^j.^: we measure a trec-ipers power 
law index in regime B of -3.2±0.5 (see further discussion 
in Sec. (H). We also calculate the mean alpha param- 
eter (a=LpersX ^rec / L'b) from the daily-averaged val- 
ues reported by ILinares et al.l ()2011D and, after applying 
the bolometric correction we obtained from broadband 
X-ray spectral fits (Sec. [U, we find a mean a of 102, 
with a standard deviation of 27. Figure |9l shows that a 
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TABLE 4 

Summary of bursting regimes in T5X2. See Section (32] for details. 



Regime 


(-^^Edd) 


^rcc 

(s) 


(10^9 erg) 


-^pcak 

(10^7 erg s-i) 


Description 


A 


0.1-0.2 


>2000 


1.6-1.3 


4.6-3.7 


Slow decrease of E^, & Lpoak 
(and tree?) with increasing Lpcra. 


B 


0.2-0.3 


2000-400 


1.3-0.3 


3.7-1.2 


Fast drop of tree, & ^poak 

with increasing Lpcrs. iroc oc Z/pcrs"^. 


C 


0.3-0.5 


400-200 


~0.3 


~1.2 


Bursts evolve into mHz QPOs with 
increasing frequency (tree ex Lpors"^). 
-Eb & ipcak "saturate". 


D 


0.3-0.5 








No bursts/mHz QPOs detected on 
normal & flaring branch in Z source phase. 



Lpci.s~0.3 L^dd (regimes A & B), and decreases at higher 
luminosities (regime C) to reach again values close to 60. 
Figure [5] also shows the burst duration and rise time as 
a function of Lpers^ the general trend is for bursts to be- 
come shorter and have slower rise when Lpcrs increases. 
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Fig. 10. — Magnetic field strength vs. spin period in accret- 
ing NSs. The dotted rectangle in the lower left corner (labeled 
'bursters') shows the location of all previously know thermonuclear 
burst sourc es with me asured spin (Psaltis iz Chakrabartv 1991; 
[Uhakrabartv et al.l 1200 3: Patruno 2010, and references therein). 
The dotted rectangle in the upper right corner (labeled 'HMXBs') 
shows the location of X-ray pulsars in high-mass X-ray binaries 
(|Bildstcn ct al.iil997i : filled gray circles show X-ray pulsars with 
_B measured from cyclotron lines, Caballerp & Wilms 2011). Red 
double-head arrows show the location of T5X2 (using the two mag- 
netic field estimates from Papitto et al. 2011; Miller ct al. 2011). 
Vertical black lines mark the location of other slow pulsars in 
LMXBs (Ps>0.01 s), none of which has shown thermonuclear 
bursts to date. These are labeled with the first two characters of the 
sourc e names, which follow: 2A 1822-371 ( Jonker & van dcr Klis 
[2001D . 4U 1626-67 (Rappap ort et al.| [T977l: Orlandini et al. 1998), 
GRO 1744-28 (the "bursting pulsar", BP; Kouvcliotou ct al. 199_^; 
[Fmger et al. 1996; Cui 1997; Sturner & Dermer 1996), Her X- 
1 (Trucmper ct al. 1978; Bildsten et al. 1997; Mihara et al. 19^; 
Idal Flume et al. 1998) and GX 1+4 ( Chakrabartv fc Roche 1991 : 
iTTmkle et al.ll200Q : lCuilll997l') . 



4. DISCUSSION 

The bursting behavior of T5X2 strikes as surpris- 
ing for several reasons. During the last three decades 



bursts at persistent luminosities L^dd bad 

proven exceptional and extremely difficult to detect, 
even when studying a large sample of bu rsters (e.g. 
iCornelisse et al.l 120031 : iGallowav et all I2008L and refer- 
ences therein). This decrease of burst rate at high 
Lpersi opposite to what standard burs t theory pre- 
dicts (iFuiimoto et al.lll981l: [Bildstenl ll99g). has been at- 
tributed in the literature to several effects, including: i) 
stable thermonuclear burning becoming more important 
at high m an d consuming an increasin g fraction of the 
accreted fuel (jvan Paradiis et aP 119881 ): ii) deflagration 
fronts or "flames" propaga ting on the N S surface and 
consuming part of the fuel ()Bildstenlll995D and iii) non- 
spherical accretion confined to a fraction of the NS sur- 
face that would increase with Lpors (Bildsten 2000) . The 
large number of bursts observed from T5X2 in October- 
November 2010 (at an average burst rate over more than 
a month of 4.9 hr~^; Sec. [3]) make T5X2 the most pro- 
lific source of thermonuclear bursts known to date. No- 
tably, such copious burst activity was observed when 
Lpors was in the range 0.1-0.5 LEdd, a regime where 
burst rate is found to decrease drastically in all other 
bursters. The only source that has shown (albeit sporadi- 
cally) such high burst rates at similar Lpcrs is Cir X-1: on 
May 2010 it featured burst recurrence tim es as short as 
1000-2 000 s, when Lpcrs was ^0.2 LEdd (Linares et alJ 
l2010bf) . Two other well known sources of "high-m" 
bursts, GX 17-1-2 and Cyg X-2, have shown bursts much 
less frequently (mean burst rate over more than 10 yr 
was 0.05 and 0.1 hr~^, respectivelv: IGallowav et aT1l2008D 
and when Lpcrs was close to Lsdd- Other systems have 
shown thermonuclear bursts at intermediate accretion 
rates (Lpors^O.1-0.5 LEdd; including GX 3-1-1, GX 134-1, 
Ser X-1, 4U 1636-53, 4U 1735-44, 4U 1746-37), but 
only sporadically ( ^wa.it, of at least one hour and typi- 
cally much longer; IGallowav et al.l [20081 and references 
therein). T5X2 is therefore exceptional in that it fol- 
lows the expected b urst rate at high m (jFuiimoto et alJ 
fimtlBiIdstml ll998l) much more closely than any other 
know n burster ([Cornelisse et alJ 120031 : IGallowav et alJ 
120081 ). We extend on this comparison between thermonu- 
clear burning theory and T5X2 in Section [43] 

Moreover, T5X2 is the first slow X-ray pulsar (spin pe- 
riod Ps>10 ms) to show thermonuclear bursts. None of 
the slower (LL^l s) "classical" X-ray pulsars in HMXBs 
have shown thermonuclear bursts to date, even if their 
persistent luminosity varies over the same range as in 
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bursters. This is usually attributed to the stronger dipo- 
lar NS magnetic field (B) in HMXBs channeling the ac- 
cretion flow into a much smaller area than LMXBs, lead- 
ing to stable bur ning of all the a ccreted fuel due to a 
very high local m (jJoss fc Lilll980l: see also Sec. l4.3p . On 
the other hand, only 7 out of the more than 90 bursters 
known have shown X-ray pulsations in the persistent 
emission (i.e., 7 out of the 14 accreting millisecond pul- 
sars have shown bursts to date) . Before the discovery of 
T5X2 the slowest spinning burster had a spin frequency 
more than 20 times higher than T5 X2 (ICR J17511- 
3057 , with a spin frequency t/s= 245 Hz: lMarkwardt et al.l 
120091 lAltamirano et al.ll2010aD . T5X2 therefore bridges 
the gap between "pulsars that don't burst a nd bursters 
that don't (typically) pulse" ()Bildstenl [19981 and refer- 
ences therein). We show this in Figure [TU] by comparing 
both the B and values of T5X2 to those of LMXBs 
(including a few peculiar LMXBs) and HMXBs. Fig- 
ure [TU] clearly shows that T5X2 features values of B and 
Ps intermediate between bursters and HMXBs. We dis- 
cuss the consequences of such a high B and lon g Pg for 
thermonuclear burst regimes in Sections 14.31 and 14. 4[ re- 
spectively. 

The shortest burst wait times known to date 
([Linares et al.ll2009t iKeek et al.ll2010l 5.4 and 3.8 min, 
respectively) were based on sets of typically two or 
three consecutive bursts followed by much longer pe- 
riods without bursts. The bursts from T5X2 pre- 
sented herein are remarkably quasi-periodic, not only 
in the mHz QPO phase, as witnessed by the well 
defined and smoothly evolving wait time in individ- 
ual bursts (twait; see Figs. [T] & [5|). This behavior 
is analogous to the "clocked burster" (GS 1826-24 ; 
Tanakal[l9 8^: Ubcr tini et all [1991 [Gallowav et al."2004'; 



Heeer et al . 2007a) and IGR J17511-3057 (Falanga et al. 



20111 ) , which have shown bursts at regular intervals with 
an approximate t^cc oc Lp^g relation, although at lower 
accretion rates than T5X2. For this reason, we compare 
in Figure [11] the burst energies and recurrence times of 
T5X2 to those of these two sources, as well as Cir X-1. 
The bursts from T5X2 (at m ~0. 1-0.5 mEdd) are more 
frequent and less energetic than those from GS 1826-24 
and IGR J17511-3057 (at m -0.02-0.06 mEdd)- Fig- 
ure [TT] also shows that the t^cc oc L~^^.^ relation that 
we find in regime B is unique to T5X2. The bursts 
from Cir X-1, however, have again a strong resemblance 
to those of T5X2, with energies in the same range 
(< 10^9 erg) and a similar i?b-ipers relation (a break 
or transition is also suggested by the Cir X-1 data; Fig- 
ure [H]). 

4.1. Bursting regimes vs. burning regimes: 
the need for heat 

In the present Section we compare in detail the T5X2 
burst properties with theoretical predictions. Models 
of nuclear burning in the envelope of (non-magnetic, 
non-rotating) accreting NSs predict four different burn- 
ing regimes on a NS accreting a mixture of hydrogen 
(H; mass fraction Xg), helium (He; mass fraction Yq) 
and heavy elements (m ainly CNO, mass fra ction Z; see 
Wooslev &: T^am|[l976l;rF uiimoto et al. 198l[; lTaam|[l98l[; 
Bildstenlll998l;lCumming fc Bildstenii2000i and references 
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Fig. 11. — Burst energy (top panel) and recurrence time (bot- 
tom panel) vs. persistent luminosity (bottom axis) and inferred 
mass accretion rate per unit surface (top axis) for four bursters, as 
indicated on the bottom panel: T5X2 (filled black circles; this 
work), GS 1826-24 ("the clocked burster", open gray squares; 
from Gallowav et al. 2004), Cir X-1 (open red pentagons; from 
[Linares et al., ,2010bi) and IG R J17511-3057 (open gray trian- 
gles: from lFalanga et aTll20lll 1. The solid line shows the empiric 
tree cc Lpe^s^ relation found for GS 1826-24 by IGalloway et al.l 
l(200l), similar to the t 

rec cx^ ^pers relation found in IGR J17511— 
3057 IjFalanga et al.12011] '). We also show the Ei,-rn and trcc-rn rela- 
tions predicted by two ignition models from Gumming & Bildste^ 
II2000I . labeled CBOO) with different accreted hydrogen fractions: 
Xq=0.7 (green dashed line) and Xo=0.1 (magenta dotted line). In 
both models the metallicity was assumed to be Z=0.02 and the 
heat flux from the crust was fixed at Q^=0.1 MeV nucleon"^. 



than mEdd: both H and He burn stably and no bursts are 
expected (Sec. 14. 2| ). For accretion rates 



m > msHb — 0.008 mEdd 



0.7 



0.01 



1/2 



(1) 



therein). At the highest accretion rates, close to or higher 



(jBildsteni 119981 where we assumed an opacity of 
0.04 cm^ g~^ and the value of rh-Edd given in Sec. [5]), H 
burns stably between bursts at a consta nt rate, via the 
so-called (/3-limited) "hot- CNO" cycle (IBildstenI [19981; 
iCumming fc Bildstenll2000l and references therein). In 
this range all bursts are triggered when He burning at 
the base of the accreted layer becomes thermally unsta- 
ble (He ignition). The column depth (or density) at the 
base of the burning layer when ignition occurs is known 
as ignition depth, yign, and the time between bursts is 
simply tree = 2/ign/m. 
As H burns at a constant rate (temperature and m 
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independent, but proportional to Z), the time it takes to 
consume all H in a sinking fluid element depends only on 
Xq and Z: 



tburn = 27 hr ( ^ 



0.01 



(l + z)/1.31 (2) 



(as seen by an observer; e.g., iGallowav fc Cummind 
Hoop), where l+z=(l-2GA//i?c2)-i/2^1.31 is the grav- 
itational redshift on the surface of a NS with mass 
M=1.4 Mq and radius i?=10 km. The longest burst 
recurrence time that we measure in T5X2 is ircc ~1 hr. 
For Solar abundances this implies that tree << iburn, 
i.e., there is no time to consume all H between the T5X2 
bursts unless the fuel is substantially H-poor and/or 
metal rich (see below). In general, for accretion rates 



0.7 

m > mdop ^ 0.04 rliEdd ( -rr 



0.01 



13/18 



(3) 



(|Bildsteni ri998. for M=1.4 Mq and i?=10 km), there is 
no time to burn the accreted H before reaching ignition 
conditions (i.e., t^oc < iburn) so that He ignites in a mix- 
ture of H and He. In this regime j/ign does not depend 
sensitively on rri, and a simple t^ec oc relation is 

expected. For m < rridcp instead, there is enough time 
to deplete all H before the base of the accreted layer 
reaches ignition conditions. Bursts are then triggered in 
the absence of H (pure He ignition). In this regime ?/ign 
decreases with increasing m, which results in a steeper 
decline of the burst recurrence time as m in creases, close 
to tree oc m~ 3 (|Cumming k Bildstenll20M ). 

This transition between pure He and mixed H/He ig- 
nition regimes can be seen in the E\,-rh and t^cc-'ni re- 
lations in Fi gure | 11[ where we sh o w tw o semi-analytic 
models from lCumming fc BildstenI (|2000) with different 
compositions. Even though the transition from pure He 
to mixed H/He ignition is expected at rridcp~0.05mEdd 
for the case of Solar abundances. Figure [Tl] and Equa- 
tion [3] clearly show that changes in the accreted composi- 
tion can increase Wdcp by a factor of at least 10. In par- 
ticular. Figure [Til shows that ignition models with low H 
abundance, [Xq—0.1, Z=0.02], can reproduce the change 
in slope of both the i?b-ipcrs and trcc-^pers relations that 
we observe in T5X2 at ipcrs — O.S-Ls dd, i -e-, at the tran- 
sition between regimes B and C (Sec. 13. 2[ Tables |3] & S]) . 
With such low Xq, H can be depleted before reaching He 
ignition at accretion rates much higher than in the case 
of Solar abundances, which highlights the importance of 
fuel composition in the burning regimes. 

The average T5X2 accretion-to-burst energy ra- 
tio (a=(5grav(l+^)/(3nuc), (a) = 102, corrcspouds 
to a total nuclear energy release during the bursts 
Qnuc=i2.8 MeV nucleon-^ (see also iMotta et al.ll201ll ). 
Taking Qnuc = 1-6 + 4(X) MeV nucleon-^ (which as- 
sumes complete burning of the accumulated fuel and 
~35% neutrino energy loss; IGallowav et al.l 120081 and 
references therein) we find an average H mass fraction 
over the burning layer (X)~0.3. Such low inferred Qnuc 
lends support to the low H fraction fuel scenario for T5X2 
proposed above. In summary, the slopes of the i?b-ipors 
and trec-Lpers relations as well as the a values observed 
during the T5X2 outburst, indicate a low accreted H frac- 
tion and strongly suggest a transition from the pure He 



ignition regime to the mixed H/He ignition regime hap- 
pening during the outburst rise when m increases above 
~0.3 rriEdd- It is worth noting that the reverse transition 
is observed during the outburst decay when m drops be- 
low 0.3 TTiEdd (i-G., the outburst rise and decay tracks 
overlap and no hysteresis is seen in Figures [7] & [H) . 

We stress that the link between the observed T5X2 
bursting regimes and the theoretical burning regimes that 
we put forward is based on the iroc-ipers and i?b-ipers 
relations. Modeling and interpretation of the individual 
burst light curves is beyond the scope of this work, yet 
we note that different burst light curves and peak lu- 
min osities are predicted i n the different ignition regimes 
(e.g. iWooslev et all 120041 ). Pure He bursts are expected 
to show faster rise times than mixed H/He bursts. This 
agrees qualitatively with the identification of regime B 
as pure He ignition: t^sc is shorter in regime B than 
in regime C (Fig. [9|). The large change in Lpcrs, how- 
ever, could influence the observed timescales along the 
burst rise and decay, as these are measured after sub- 
tracting the persistent emission. It should also be noted 
that at low accretion rates (near TOsHb — 0.08 WEdd 
for Xq=0.1, Z=0.02) there should be a transition to H- 
ignited burst s, not included in the mod els discussed in 
this Section (jCumming fc BildstenlfeOOOD . 

If the link between bursting and burning regimes that 
we propose is correct, to our knowledge this is the first 
time that the transition between pure He and mixed 
H/He ignition is observed in a single source. There 
is, however, a systematic and interesting discrepancy 
evident in Figure 1111 Even when including compres- 
sional heating or when increasing the base heat flux 
from the NS crust (up to 2 Me V nucleon"^; fixed at 
0.1 MeV nucleon" -1 in Fig. [Til [Gumming fc BildstenI 
|2000[) . ignition models predict higher burst energies and 
longer recurrence times than those we find in T5X2, by 
a factor close to 10 in most cases (i.e., larger than the 
distance uncertainty on Lpors and i5b, Sec. [2]). This large 
difference between the observed and predicted values of 
£'b and tree, together with the lack of detailed model- 
ing of the T5X2 burst light curves, prevents a conclusive 
identification of the observed bursting regimes. 

We propose that such discrepancy could be explained 
by the presence of an extra source of heat in the NS 
envelope not accounted for by ignition models (which 
typically consider only hot-CNO heating). Additional 
heat would act to reduce yign and thereby decrease E'b 
and iroc, explaining the T5X2 observations presented 
herein. Several interesting possibilities for the nature 
of this extra source of heat have been partially inves- 
tigated, including: i) heating due to stable He burn- 
ing (triple a lpha reaction rate s have been recently de- 
bated; e.g., 
IPeng fc OttI 



Ogata et al.l [20091: [Potter fc PaxtonI [20091: 
20101): ii) heating du e to deep burning of 



residual H (e.g.. lTaam et al.lll996L who already pointed 
out that it can lead to weaker and more frequent bursts 
than expected); iii) thermal inertia, "hot ashes" or heat 
from previous bursts, which could become important for 
T5X2 as it features ircc ~ 200s, the shortest recurrence 
times between thermonuclear bursts ever observed (time- 
dependent simulations of a series of bursts are needed 
to investigate this in detail; e.g., iHeger et al.l l2007b[ ) 
and iv) turbulent friction at the base of the spreading 
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layer, which could release s ubstantial amounts of heat 
(|Inogamov fc Sunvaevll2010l ). It is wort h noting that an 
indep endent study of the same source ()Degenaar et alj 
1201 It ) reached similar conclusions, suggesting the pres- 
ence of a "shallow heat source" in T5X2 in order to rec- 
oncile quiescent observations and crust heating/cooling 
theory. Without regard to the exact nature of this heat 
source, we have shown in the present work that the burst 
properties of T5X2 place new constraints on the thermal 
properties of a fast- accreting and frequently-bursting NS. 

4.2. mHz QPOs and marginally stable burning 

The mHz QPOs from T5X2 (Sec. O see also 
iLinares et al.l l2010af) have distinctive properties that 
clearly set them apart from the previously know n mHz 
QPO s (Sec. [TJ Rcvnivtscv et al. 2001; Alt amir an o et al.l 
120081 ). The persistent luminosity that we measure in 
T5X2 when mHz QPOs are present is about 10 times 
higher than that observed in previous mHz QPO sources 
(Lpors higher by a factor of 4-25 taking into account 
the observed ranges: 0.02-0.1 L-^dd hi atoll sources as 
opposed to 0.4-0.5 iEdd in T5X2; Table [2]). For this 
reason we refer to the previously known mHz QPOs as 
"low-Lpers mHz QPOs". 

Bright bursts and low-Lpors mHz QPOs alternate, 
while remaining clearly distinguishable. Instead, in 
T5X2 bursts smoothly evolve into mHz QPOs and vice 
versa (Sec. 13. 1|) . a phenomenon never observed before. 
Strikingly, the same qualitative evolution from bright in- 
frequent bursts to faint and frequent bursts, mHz QPOs 
and ultimately stable burning is predicted to happen 
as m increases by both one-zone models and detailed 
simulations of nuclear burning on NSs accreting near 
the boundary betw een unstable and stable He burning 
(|Heger et aim mh. We can therefore identify with con- 
fidence the mHz QPOs from T5X2 with marginally stable 
burning on the NS surface. The evidence that links low- 
Lpers mHz QPOs with marginall y stable burning is less 
conclusive, but remains valid f sec Rcvnivtscv ct al.l l200ll : 
lYu fc van d"erKlis 2002: Ahamirano ct al. 2008, for de- 
tails). 

Analytic estimates place the threshold for stable He 
burning at 



1.7 



1 + Xr 



3/4 



0.3 X 0.6 



1/2 



(4) 



wher e At is the mea n molecular weight of the accreted 
fuel (IBildstenI [l998L assuming again M—IA Mq and 
i?=10 km). 

The persistent luminosity where mHz QPOs are ob- 
served in T5X2 (0.4-0.5 Lsdd) is therefore closer to the 
expected value of rrisb than what was seen in low-Lpds 
mHz QPOs, yet still inconsistent with the value pre- 
dicted by theory, which is higher by a factor of ~2 (unless 
Terzan 5 is at ~9 kpc instead of 6.3kpc, Sec. [2]). Invok- 
ing H-poor fuel makes this discrepancy even larger, as it 
can increase the expected rhsh by a factor 2 (for Xq—0.1, 
Z=0.01; Eq. 21). A more massive NS can have unstable 
burning at higher m than a less massive star, but given 
the weak (M^/^) scaling of rhsh different NS masses can- 
not explain the large difference (factor 4-25) in Lpers be- 
tween T5X2 and the other known sources of mHz QPOs. 
Another remarkable difference betwen the mHz QPOs 



presented herein and the low-Lpors mHz QPOs resides 
in their fractional rms spectrum (Sec. [2]): we find a 
fractional rms amplitude that increase s with energy be- 
tween 2 and 20 keV (Sec. 13. 1|) . while iRevnivtsev et al.l 
(2001*) reported fractional amplitudes decreasing with 
increasing photon energy between 2 and 5 keV. The 
energy-averaged fractional amplitudes are, however, sim- 
ilar (1.3-2.2 % in T5X2, this work; 0.7- 1.9 % in low- 
Lpcrs mHz QPOs. IRevnivtsev et al.]|2001[ ). The increase 
of fractional amplitude with energy that we find is shal- 
lower than the linear increase predicted by models of tem- 
perature oscillations at the NS surface (developed in the 
context of burst oscillations: iPiro fc Bildstenll2006l ). 

An interesting possibility is that the low-Lpers mHz 
QPOs trace the boundary of stable H burning (Eq. [ij, 
whereas the mHz QPOs that we discovered in T5X2, at 
a higher accretion rate, occur at the boundary of stable 
He burning (Eq. \^ . This remains speculative at present 
given the lack of published theoretical work on this par- 
ticular topic. A specific analysis of oscillatory burning 
at the thermal stability boundary of H burning is needed 
to address this hypothesis. The oscilla tory be havior at 
the marginally stable point is generic (|Paczvn ski 1981) 
and it could also occur when H burning stabilizes. If 
low-Lpcrs mHz QPOs do happen at the H burning sta- 
bility boundary, the need to invoke confined accretion to 
e xplain their low 



would vanish. 



iHeger et all (|2007br ) found that the QPO frequency is 
mainly sensitive to the accreted H fraction (Xq) and 
the NS surface gravity: increasing the surface gravity 
or decreasing Xq leads to higher mHz QPO frequen- 
cies. The highest z/qpo values that we find in T5X2 
(4.2 mHz) are about a factor of 3 lower than those of the 
low-X n^rs mHz QPOs (Table E] IRevnivtsev et al.l 120011 : 
Altamirano et al. 2008). If the accreted fuel has a sim- 
ilar composition, and if the mHz QPOs have the same 
origin, this would suggest a less compact NS in T5X2 
than in the low-Lpers mHz QPO sources (the "atoll" 
sources 4U 1636-536, 4U 1608-52 and Aql X-1). Compar- 
ing to models of marginally stable burning, we find that 
if Xo<0.7 the surface gravity of the NS in T 5X2 must be 
q = G A//i?^<1.9xlO^'^ cm s'^ (Figure 9 in IHeger et all 



l2007b[ ). Even though it is subject to theoretical (based 
on analytic one-zone model) and observational (a new 
outburst of T5X2 could show higher z^qpo) uncertainties, 
this illustrates how a comparison between marginally sta- 
ble burning models and mHz QPO properties can be used 
to constrain the NS compactness. 

Despite the striking similarities between the T5X2 
burst properties and the general bursting behavior pre- 
dicted by models of nuclear burning near the transi- 
tion from unstable to st able burning (cf. Fig. [T] in 
this work and Fig. 5 in ' Heger et al.l l2007bl ). interest- 
ing differences remain. First, as explained above, the rh 
where marginally s table burning is expected (0.925 msdd 
IHeger et al.ll2007b() is higher than the highest inferred m 
wher e we observe mHz Q POs in T5X2 (0.5 mEdd)- Sec- 
ond, ' Heger et al.l (|2007U ) find a sharp transition from 
bursts to mHz QPOs and finally stable burning (oc- 
curring between 0.923-0.95 mEdd), while we observe in 
T5X2 a smooth evolution from bursts to mHz QPOs, 
and vice versa (between 0.1-0.5 rriEdd; Sec. Time- 
dependent simulations of nuclear burning tailored to 
T5X2 would be of high interest, in particular exploring 
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the H-poor fuel range. Third, we find that the disap- 
pearance of the mHz QPOs in T5X2 is not linked to an 
increase in to but to an actual drop of Lpcrs that hap- 
pens when the accretion state changes (from horizontal 
to normal & flaring branches). This suggests that the 
geometry /configuration of the accretion flow plays a role 
in setting the nuclear burning stability boundary. 

4.3. The role of an intermediate magnetic field 

While most models of thermonuclear bursts on ac- 
creting NSs assu me that the NS rnagnet ic field is 



negligible ( Fui imo to et al.l 
fl98l Imoslev et al.ll2004lv 



19811: iTaamlfTop iPaczvnskil 
Joss fc Lil (|1980[ ) showed that 
the presence of a strong {B>10^^ G) magnetic field 
can affect the stability of nuclear burning in different 
ways. A strong magnetic field reduces the (conduc- 
tive and radiative) opacities, allowing for more efficient 
heat transport and thereby stabilizing burning. Ap- 
plying disk-ni a gneto sphere interaction models to T5X2, 
liPapitto et all ()2011h estimated 5=2x10^-2.4x101° G 
from the lumin osity range at which pulsations were de- 
tected, whereas I Miller et al.l (|201lD further constrained 
B=(0.7-4)xl09 G from the Fe line profile. The fact 
that T5X2 shows thermonuclear bursts suggests that the 
magnetic field needed to suppress the thermal instabil- 
ity and quench thermonuclear bursts must be greater 
than ^10^° G, in acc ordance with theoretical expecta- 
tions (Jos s fc LillT980 [). 

A strong magnetic field also affects cony ec tive heat 
transport and mixing (e.g. iJoss fc LH llQSOD . iBildstenI 
(1995i) proposed that even in cases where thermonuclear 
burning is unstable the presence of a strong magnetic 
field can suppress convection in the NS envelope, stalling 
the propagation of the burning front (diffusive heat trans- 
port being slower than convection) and preventing the 
fast ignition t hat causes type I X-ray bursts (see also 
IBildsten|[T998l) . A gain, the mere presence of type I X-ray 
bursts in T5X2 strongly suggests that the magnetic field 
required to suppress convective burning fronts must be 
higher t han ^10^° G, in agreement with analytical es- 
timates (|Joss k Lil 1980; Bildstcn & Brown 1997). Rise 
times are rather long though (Fig. [H]), which could per- 
haps indicate t hat B is strong enough to slow down the 
burning fronts ()Bildsten fc Brownlll997D . In summary, 
the bursting properties of T5X2 presented herein (in par- 
ticular combined with refined measurements of its mag- 
netic field strength) can place new constraints on the 
physics of convection and heat transport under interme- 
diate magnetic fields (lO^-lQi" G). 

Finally, a 10^-10"'^° G magnetic field can also chan- 
nel the accretion flow onto the magnetic polar caps, as 
suggested b y theoretical mod els of disk-magnetosphere 
interaction (jLamb et al.lll973D and by the presence of X- 
ray pulsations in T 5X2 (jStrohmaver fc Markwardtll2010l : 
iPapitto et al.ll201ll ). Accretion confined to the magnetic 
polar caps may break the spherical symmetry assumed 
by most thermonuclear burst models, and increase to for 
a given (observed) Lpas (Sec. [5]). The impact on the 
burning properties depends on the size of the polar caps 
and on how soon the accreted fuel spreads laterally while 
sinking into th e NS envelope under the influence of a 
magnetic field. iBildsten fc Brownl (jl997| ) found that for 
dipolar magnetic fields B< (2-4) xlO^" G, the accreted 



fuel spreads before igniting and the spherically symmet- 
ric case is recovered. This suggests that the T5X2 burst 
properties can still be compared to spherically symmetric 
ignition models (Sec. 14. 1| ). 

As the size and "depth" of the polar caps are nev- 
ertheless ill-constrained quantities, we explore a simple 
scenario in which the fuel in T5X2 is confined into a 10% 
of the NS surface down to ignition depth, yign- From 
geometrical considerations, assuming the same y^^^ pre - 
dicted by ignition models (jCumming fc BildstenlHoOOl ). 
in this confined scenario bursts would be 10 times less 
energetic and 10 times more frequent than in the spheri- 
cally symmetric case, which would explain the mismatch 
between predicted and observed and tree (Sec. 14. ip . 
However, the accretion rate per unit area would also in- 
crease by a factor 10, which for the observed ipcis implies 
that unstable burning would operate at m as high as ^5 
times TOEdd- This is well above the predicted thresh- 
old for s table steady b urning of H and He (rfigb , Eq. SJ 
see also IBildstenI |1998() , and we therefore consider this 
scenario unlikely. If the accreted material contains no 
H, however, TOgb can be several times larger than TOEdd 
(|Bildsten. .l998). 

4.4. Rotation, mixing and burning 

The unique bursting behavior in T5X2 is, in essence, 
much closer to what theory predicts than any other 
burster known to date. Theoretical burning regimes were 
mostly based on models that assume a non-rotating NS 
and r adial infall of the accreted mass ([Fuiimoto et al.l 
[19811: IBildstCTl[T998h . These are admittedly simplistic 
assumptions, as NSs in LMXBs can have spin frequen- 
cies in excess of 600 Hz and the orbital frequencies in 
the innermost accretion disk are well above 1000 Hz. 
The accretion of matter with angular momentum onto 
a spinning NS introduces a shear in it s surface layers 
(jPiro fc Bildstenll2007t [Keek et al.l 120091 and references 
therein) . The main effect of rotation which may have an 
impact on burning regimes is turbulent mixing. In par- 
ticular. Keek et al. (2009) found that such rotationally 
induced turbulent mixing stabilizes He burning, decreas- 
ing the burning stability boundary (TOsb)- 

Motivated by the stronger magnetic field (factor ^10) 
and slower spin (factor >20) of T5X2 compared to the 
rest of bursters, we speculate that rotationally induced 
turbulent mixing is what sets T5X2 apart from the rest 
of thermonuclear burst sources. The stronger B field 
could channel a large fraction of the accreted matter 
into the NS magnetic poles and produce a more radial 
infiow which would then spread over the whole surface, 
minimizing shear and turbulent mixing. Therefore 
T5X2 seems to meet to a greater extent the assumptions 
of zero spin and radial infall mentioned above, which 
can explain the better agreeme nt with theory. This in 
turn suggests tha,t, as arg ued bv lPiro fc BildstenI pOOTi ) 
and iKeek et al.l ()2009( ). the effects of rotation must be 
considered to explain the behavior of most bursters. 



5. SUMMARY AND CONCLUSIONS 

We have presented the discovery of mHz QPOs from 
the NS-LMXB and 11 Hz X-ray pulsar T5X2, as weU 
as the full bursting properties during its 2010 outburst, 
when the persistent luminosity varied by about a factor 
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of 5 (between Lpcrs— 0.1 0.5 LEdd)- T5X2 showed copi- 
ous thermonuclear bursts when Lpors^O.2 £Eddj a regime 
where thermonuclear bursts had proven exceptional to 
date. Burst energies and recurrence times gradually de- 
creased as I/pers increased, turning into a rapid series 
of faint bursts and smoothly evolving into a mHz QPO 
near the outburst peak. Only at the highest Lpcrs range 
bursts became undetectable during a few short intervals. 
This behavior is both unprecedented among bursters and 
remarkably similar to the nuclear burning regimes ex- 
pected on an accreting NS near the boundary of stable 
He burning. 

We find foiir different bursting regimes when studying 
the relation between burst properties and Lpers and show 
that in one of such regimes the burst recurrence time de- 
cays steeply with increasing ipers, close to tree oc L~^^g. 
By confronting the change in burst properties (in partic- 
ular, the relation between i?b, iroc and Lpcrs) to ignition 
model predictions, we find evidence of a transition be- 
tween pure He and mixed H/He ignition occurring in 
T5X2 when Lpers— 0.3 Lsdd- We note that large discrep- 
ancies remain between the observed T5X2 burst proper- 



ties and those predicted by theory at high mass accre- 
tion rates. We further argue that the accreted fuel is 
H-poor and suggest that an additional source of heat in 
the NS envelope is needed to reconcile the observed and 
predicted burst properties. 

We examine the properties of the mHz QPOs from 
T5X2 in the context of marginally stable burning mod- 
els, and compare them to those of previously known mHz 
QPO sources. T5X2 features mHz QPOs with lower 
frequencies (by a factor ~3), when Lpers is substantially 
higher (by a factor ~4-25). Finally, we discuss the role of 
magnetic field and spin in setting the unique T5X2 burst 
and mHz QPO behavior, and speculate that the absence 
of rotation effects such as turbulent mixing of the ac- 
creted fuel may set T5X2 apart from the rest of bursters. 

We thank M. van der Klis for detailed comments on the 
manuscript. We are grateful to the International Space 

Science Institute in Bern, where part of this work was 
completed. ML acknowledges support from the NWO 
Rubicon fellowship. 
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